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Abstract 

The Mean Square Displacement is a central tool in the analysis of Single Particle Tracking experiments, 
shedding light on various biophysical phenomena. However, as we show, it suffers from two systematic 
errors when analysing tracks of anomalous diffusing particles. The first is significant at short time dif- 
ferences and is induced by measurement errors. The second arises from the natural heterogeneity in 
biophysical systems. We show how to estimate and correct these two errors and improve the estimation 
of the anomalous parameters for the whole particle distribution. As a consequence we manage to char- 
acterise ensembles of heterogeneous particles even at very short and noisy measurements where regular 
time averaged mean square displacement analysis fails. This procedure has the potential to improve 
experimental accuracy while maintaining lower experimental costs and complexity. 

Notation 

Anomalous diffusion constant - Da 

Apparent anomalous exponent due to heterogeneity - as 
Apparent anomalous exponent due to noise - ajv 
Apparent location - 

Apparent TAMSD - J| 

Control experiment noise MSD - Nc 

Dynamic functional - tp 

Error in anomalous exponent due to heterogeneity - da 
Mean Logarithmic Square Displacement (MLSD) - A(a) 
Mean of anomalous exponent distribution - ^a 
Noise and heterogeneity corrected MSD - i^ciA) 
Noise Corrected MSD (NC-MSD) - i/(a) 
Noise MSD - N 

Relative error caused by noise - eat 

Standard deviation of anomalous exponent distribution - ctq. 
Standard deviation of noise - p 
Time averaged MSD - Sj^ 
Time difference - A 
True location - 

True particle anomalous exponent - a 



Introduction 

Single Particle Tracking (SPT) is an increasingly popular technique in the study of biophysical systems [l]. 
Typically, a single particle is tracked as it moves and its trajectory is extracted. This measurement can 
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be performed both in vivo or in vitro. Advances in biochemical labeUing, nano particle production and 
the use of fluorescent proteins now allows to label almost any biological entity such as, single stranded 
RNA in bacteria j2j, nano-beads in a polymer solution |3j, membrane dyanmics [4], and specific nuclear 
entities (5j|6]. 

At each time point, the location of the particle is measured and stored for later analysis. In most 
biophysical systems, the motion of the tracked particle is stochastic - thus further analysis is necessary 
before claims regarding the underlying nature of the system can be made. 

Perhaps the most basic analysis technique is the Mean Square Displacement (MSD) in which one 
calculates the average square of the distance x{t) past by the particle in time t. This can be performed in 
two ways. The first is the ensemble averaged MSD (EAMSD), in which the averaging is performed over 
the displacements of N particles, at some t since the beginning of the experiment, < >(t)~ ■ 
The second is the time averaged MSD (TAMSD) in which all displacements over a time difference of A 
from a single trajectory are averaged over: 
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Stochastic processes can be characterized according to their MSD behavior, be it the EAMSD or the 
TAMSD. If the MSD is linear in time, i.e < >= Dt, then the process is termed normal diffusion, 
where D is the diffusion constant. This is the behavior of, for example, Brownian motion. Other MSD 
behaviors are termed anomalous diffusion 't'.'s]. 

If the MSD is described by a power law, < x'^ >— Dat" then a further distinction is made. When 
1 > a > the process is called Subdiffusion and for a > 1 the process is called super diffusion. Da 
is the anomalous diffusion constant and a is the anomalous exponent. We stress that the TAMSD and 
the EAMSD do not have to show the same behavior, and that the behavior of a process may change 
depending on the analysis time t (or time difference A) j9|. 

The physical and mathematical origins of anomalous diffusion are diverse. The most popular models 
are the Generalized Langevin equation (GLE) and Fractional Brownian motion (FBM) which can describe 
the motion of many body systems and viscoelastic media 10 ■ 14 ; Continuous Time Random Walks 
(CTRW) which can describe biological trapping [9 ,15 16 ; and fractal landscapes which describe, for 



example, highly obstructed motion 17 18 



The MSD itself cannot be used to distinguish between these models alone as they can all lead to 
similar anomalous exponents. This raises a significant challenge when analysing biophysical data. Since 
in many systems SPT is applied in order to probe the dynamics and physical mechanisms governing the 
system, the ambiguity of the MSD results hinders clear conclusions. However, with the advancement of 
mathematical research, more subtle tests have been proposed which succeed in identifying the nature of 
the stochastic process behind the random motion 19]- 23 . Some of these tests have been implemented 
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successfully to experimental and simulated systems [4 13 

Another important step when studying the biophysical system causing the stochastic motion is the 
extraction of accurate parameters that characterize the system. Accurate parameters can shed light on 
the physical process taking place. For example, the dynamics of various polymer models can be described 
by a GLE model, but each will give rise to a different anomalous exponent. As an example, a phantom 
Rouse chain will have a — 0.5 while a Zimm model gives rise to a = 2/3. Various chain properties can 
alter this value 26 , and in the long time domain, all chain models behave similarly according to the 
crowding of the media 27 . Thus, by identifying both the mathematical model and the correct anomalous 



parameter, one can gain deeper understanding of the physical process. 

This task is not a simple one. Biological systems are characterized by heterogeneity of both the 
tracked particles and the surrounding medium. In addition, the total time that the biological sustem can 
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be measured is limited due to photo-bleaching, induced photo-damage to the system and cell migration. 
This inserts a variation of the motion parameters between the various particles, even in the same system. 

Since the system is heterogeneous, the average tracer dynamics and the spread of parameters is 
important. If the single particle TAMSD is noisy, and cannot be analysed accurately than the ensemble 
averaged TAMSD, ^'^a^i must be calculated and fitted to its parameters, hopefully reproducing the 
average behaviour. In this averaging, many more data points are available for each time difference than 
by taking the time average of each particle alone. Evidently, random errors are significantly reduced. 
However, one cannot accurately asses the distribution of parameters with this technique. The assessment 
of this distribution is essentially left dependant on the accuracy of individual TAMSDs. 

This challenge in the assessment of individual parameters has been found regarding the diffusion 



constant in normal diffusion 28 -30 . In these studies, the distribution of fitted diffusion constants and 
their dependence on measurement rate, track length and more have been developed. However, to our 
knowledge, no such framework has been developed for anomalous diffusion. Specifically, an understanding 
of the systematic errors of the above techniques is missing. 

In what follows we take an experimentally applicative approach regarding the accurate extraction of 
the distribution of anomalous exponents from experimental data. We identify two systematic errors that 
are expected to arise in many biophysical SPT measurements. The first is an offset due to the limited 
spatial precision (location accuracy) and the second is due to the heterogeneity of the individual particle's 
behaviour. 

As we show, both can be estimated much better by using common mathematical methods. We intro- 
duce two improved measures that substitute the standard MSDs: 1) the Noise Corrected MSD (NCMSD) 
and 2) the Mean Logarithmic Square Displacement (MLSD). With the combination of these two measures, 
the average and standard distribution of anomalous exponents can be estimated correctly. In addition, 
the average particle TAMSD can be calculated to give the average diffusion constant. Simulations of 
anomalous motion illustrate both the systematic errors that occur while using simple MSD averaging and 
the implementation of the improved method. 

With these simulations we show how to characterise very short and noisy trajectories, an impossible 
task with regular MSD analysis. As we explain in the conclusions, this new capability has the potential 
of reducing experimental costs and complexity while at the same time improving the accuracy of the 
results. 



Discussion and Results 
Short time a deviation 

Every measurement has inherent noise and uncertainty. In SPT measurements, the location X[f) can 
only be measured with a varying noise term iV(j) that depends on various factors including the point 
spread function (PSF) of the optical system, the acquisition rate, and other optical and mechanical 



limitations 30 31 . In general, N^i-^ is a random variable and its distribution dependent on t and x^j). 
However, we consider the simpler case for which the noise is not dependent on time, i.e Pjqify = Pn- This 
is the common situation for good SPT experimental conditions - negligible tracer bleach and constant 
spatial optical limitations. 

We further assume, based on the central limit theorem, that P{N) is a normal distribution with zero 
mean and p standard deviation. Thus each measurement is actually the sum of the true position and the 
normally distributed noise: X(^t) — X{t) + ^(t)- The measured TAMSD is therefore: 

n-A/r 

- — t^^i 
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Where n is number of time points along the trajectory, A is the size of the sHding window for averaging 
and T is the acquisition rate. 

The true process and the noise are uncorrelated. In the case where A is small and there are many 
averaged points, the cross terms between the two differences is eliminated. Also, the noise is independent 
of time and its MSD is therefore equal io N — . Thus we are left with 5\ = 5\-\- N, where is the 
MSD of the true process. For an anomalous process with diffusion constant Da and anomalous exponent 
a this gives 5'^ = D^A" + N. 

In order to extract the anomalous exponent one calculates what is commonly termed the dynamic 
functional, ip, by taking the derivative of the logarithm of the MSD according to the logarithm of time, 
d log (5^ 

(fi — — — . For sub- and super-diffusion, one expects ip ~ a. However, in the case described here, the 

olog A 

added noise term will give an erroneous anomalous exponent ,ip — a^: 

aN(A)^a (3) 

l + iV/(Z3A") 

The relative error of the measured dynamic functional is e^r = (ajsi — a) /a. After a short manipulation 
we find: 

ejv = ^ ^ (4) 

l + i6l~ N)/N 

From equation|4]and |3]several facts emerge. First of all, the noise always works to lower the anomalous 
exponent. This is in line with previous findings which have shown a lowering of the anomalous exponent 
in simulations [32| . Second, the error is time dependent. At very short time intervals, A — >■ 0, we find 
Eat — > — 1 and aN(A) ^ 0. At long averaging intervals, the error becomes negligible and aN ^ a. 

Figure ^ shows the TAMSD of a single particle with a = 0.5 and with noise of magnitude p = 1.2. 
The deviation from the theoretical prediction is seen clearly at small A values. Notice that as this is a 
single particle track, it is not smooth. The prediction of equation |3] is presented in figure [Tj3 together 
with the average dynamic functional calculated from ten individual tracks and a clear fit is seen. We 
stress that averaging a large quantity of tracks will not eliminate this systematic error as it does not arise 
due to limited sampling. 

Luckily, this systematic error can be approximately corrected with simple experimental measures. 
Suppose a control measurement is performed in order to asses the magnitude of the noise. This is 
actually a necessary step in any experimental system and can be performed by various means. In live 
cell measurements, for example, the measurement noise is sometimes classified by measuring a chemically 
fixated cell [s] . The static noise data can then be analysed to give the MSD of the noise alone -Nc which 
should be constant in time. Note that this measurement does not estimate effects such as motion blur. 



However, these are typically much smaller in SPT experiments, especially with short exposure times 30 
After finding the MSD of the noise, one can calculate the Noise Corrected MSD (NC-MSD), i^(a) = 

(5^ — Nc- Since the control experiment is assumed to give the same noise as in the actual measurement, 
one finds i^(a) ~ DA". Now, calculation of the dynamic functional will give a with only negligible error. 

This correction is applied to the previous simulations and presented in figure [TJV and B. It is clear that 
the NC-MSD accurately retrieves the anomalous exponent of the simulated process even when the noise 
is of the same magnitude as the process MSD. This experimental scenario was previously impossible to 
study without large errors. From the theoretical point of view, the NC-MSD enables the measurement 
and analysis of anomalous diffusion down to unlimited short measure rates. 
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Large time a variation 

When calculating the TAMSD of a single particle, for large A values or short tracks, there are not many 



points to average over. Hence the TAMSD suffers from high variation and random errors 28 . These 
errors are obvious for example in Figure [2]A where it is impossible to analyse the TAMSD after the first 
few time points. 

A possible way to mitigate this problem in the case of normal diffusion is to take the ensemble average 
of many TAMSDs (EA- TAMSD) of different tracks, {^"ji^- This method works adequately for normal 
diffusion or for identical particles, retrieving accurate average diffusion parameters for the system. 

However, for the general case of anomalous diffusion, where the tracers and their environment are 
not identical this method fails. In other words, the ensemble average TAMSD does not represent the 
behaviour of the typical particle. As we show shortly, the dynamic functional extracted from ^<^a^ 
different from the average anomalous exponent and shows a time dependency. 

Let us begin with the simple case of two particles. We give both of them a diffusion constant of 
D]^ = 1 and assume ai > a2- It is clear that the average anomalous exponent is "^^"^ . However, the 
average TAMSD is ^ ^"^^ ^ . Upon extracting the dynamic functional from this average TAMSD, one 
will find that at long times, i.e. A >> 1 the average MSD behaves according to ai while at A << 1 
the behaviour is governed by (/? « a2 • 

We now analyse a large ensemble of particles and take the common experimental situation where 
for each j'th particle there is a different underlying anomalous exponent, ai ~ P{oi). For clarity and 
simplicity, we take P to be a distribution with mean /i^ and standard deviation a a- We note that since 
ai must be positive, such a distribution must be limited at zero, i.e. ai > 0. However, this limitation 
is negligible unless a a — fJ-a- The relevant corrections of the following equations for this scenario are 
presented in the SI. We also take Di ~ -P(^) with a an average D. However, D^^i and ai are not 
correlated. As we will see later, the diffusion constants do not matter in the case of accurately extracting 
the anomalous exponent. This stems from our interest in the power law behaviour in long and short 
times, dimming the diffusion constants influence negligible. 

For the above distribution of anomalous exponents, with each trajectory behaving according to the 
theoretical prediction 5^ — 2Z?ct ,jA"% then the ensemble average TAMSD behaves according to: 

'5l)^DA^--exp[ialHAf)/2] (5) 

It is clear that there is a time dependent offset from the average anomalous exponent /Iq . In order to 
find the time dependent anomalous exponent, ^^(A), we again calculate the dynamic functional {(p — as). 
One finds an effective average anomalous exponent of: 

as = + crl log(A) (6) 

From equation |6]we find the difference in the dynamic functional between two time lags. For Ai and 
Ao, and defining da = a(Ai) — o:(Ao) receive: 

da = allog{^) (7) 

Notice that da is independent of the diffusion constant and shows a logarithmic dependency on Ai. 

To test this result, 1000 particles with a distribution of individual anomalous exponents according to 
/ic = 0.5 and ct^ = 0.2 where simulated for 1000 time points. The retrieved distribution of anomalous 
exponents is shown in Figure [2j3. Two graphs of the MSDs of individual particles are shown in Figure 
[2}\. The EA-TAMSD of the ensemble is shown in Figure [2p. Notice that the time difference and MSD 
are plotted on a logarithmic scales and a sub-diffusive MSD is presented as a straight line with slope a. 
Evidently, the EA-TAMSD has a changing slope, as expected from equation |6] The dynamic functional 
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extracted at each time point is seen in Figure [2jl). It is clear that the error in the dynamic functional 
grows linearly with the logarithm of the time differences. Notice that these simulations where preformed 
without the addition of measurement noise, and thus the single source of systematic error is the spread 
of anomalous exponents. 

The time dependency of the dynamic functional is the attempt to compare between different powers of 
time. Hence, one should try to directly average the exponent values while cancelling the exponential time 
dependency of each trajectory. This can be done by averaging between the logarithm of the exponential 
quantities - creating a linear dependency in time. For this purpose we present the Mean Logarithmic 
Square Displacement (MLSD): 

C(A)-iog[E (8) 

For a process that shows anomalous TAMSDs, The MLSD gives for each particle i, ^(a) = log[£'i] + 
ai log[A]. Now, when taking the ensemble average, the true average exponent is found - 

(C(A)) = (log(i?))+Malog(A) (9) 

And specifically, (p = /ia for all times. The results of the MLSD for The same ensamble of particles in 
Figure [2J3 is also presented in Figure [2p and D. The perfect reconstruction of a™ can be seen. 

Extracting the standard deviation of 

The MLSD gives more details about the anomalous behaviour of the diffusing particles than just ^a- 
One can also extract the standard deviation of the individual anomalous exponents, (Tq.. 

At each time point, where the slope of the MLSD is taken q;(A), one also takes the slope of the 
regular MSD, as{A). Then, extract a set of the errors — a{Ai) — a5(A,j) in estimation of the average 
anomalous exponent at times A^. This set can then be fitted against Log(Ai) according to a small 
manipulation of equation [7] — cr^ log(Ai) — cr^ log(Ao). The slope of this fit, is the variation of the 
distribution of anomalous exponents in the population of particles, cr^. 

It is important to notice two central points regarding this technique. The first is that both the MLSD 
and MSD are effected similarly by the noise induced error, e^r. Thus the A^'s used for the fit can extend 
to small A values, even though might be significant. This greatly increases the amount of points that 
can used for the fit. 

Second, in some experimental scenarios, this technique can be much more accurate than fitting each 
individual TAMSD to a power law. For example, if only a small amount of time points are measured, 
the individual TAMSD is very noisy. Thus the ability to extract individual anomalous exponents from 
each track is very low. 

In order to test this technique, 5000 particles with ai taken from a normal distribution with mean 
fia = 0.6 and standard deviation a a = 0.15 where simulated. To further complicate the analysis, a 
measurement noise of p = 1.5 was incorporated into the process. For each track, only 2^ time points 
were taken, thus the ability to fit each anomalous exponent is very low. Figure [3]A_ shows the histogram 
of the simulated values and the fitted ones. Two errors are evident. The first is the very wide spread 
of the fitted distribution, and the impossible negative values found in some cases. Actually, the standard 
deviation of the fitted a^'s was cr^** = 0.36, an error of over 100%. The second error is the shifting of the 
peak value do to the measurement noise to an average value of 0.4. 

Figure [3J3 presents the extraction of CTq with the use of the MSD and MLSD. The measurement noise 
causes the dynamic functional found for the MLSD to be well below the average value of /io, — 0.6. 
This however does not effect the result, as the standard deviation of alpha values was found to be 
^MLSD _ Q jg^ g^j^ error of only 7%. We see that the MLSD-MSD difference technique extracts CTq, even 
for short, noisy trajectories, a task not possible previously. 
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Combining both techniques 

With the tools of the NC-MSD and the MLSD, one may construct the average particle MSD from 
the measured trajectories. In doing so, one must take into account that the NC-MSD is affected by 
the distribution of anomalous exponents and will suffer an error according to equation [s] dimming it 
inadequate for large A. Since the MLSD cannot be easily corrected for measurement noise (as the 
logarithm cannot be negative) it will suffer an error according to equation |4] and give < Log{Di) > 
instead of < A >. 

The key is to extract CTq from the distribution induced error with the use of the MLSD-MSD differences 
as described above. Then, since we obtain an estimation for equation [Tj the NC-MSD can be corrected 
to the true average dynamic functional by division. 

^c(A) = ^ (10) 

This procedure was implemented for a set of 1000 particles, measured at 1024 time points, with 
fia = 0.6 and — 0.2. In addition, a measurement noise of p = 1.5 was added. Here too, the analysis 
of individual trajectories is bound for random errors. 

First, the ensemble distribution error, da was extracted by comparing the MSD to the MLSD as 
stated above. Then, the NC-MSD, i^(A), was calculated for each time point. Finally, equation 10 was 

used to find i'c(A). The uncorrected EA-TAMSD, (^^a^ ^^c(A) are both presented in Figure 

It 

with 

the theoretical average particle MSD. The significant improvement is clear. Finally, Figure [4|3 shows the 
dynamic functional at each time point of both the uncorrected MSD and the final noise and variation 
corrected MSD. Evidently, equation[lO]reproduces the average particle MSD accurately - both the average 
diffusion coefficient and average anomalous exponent are reproduced. 



Conclusions 

In this work we have studied two types of systematic errors in the estimation of anomalous diffusion 
MSDs, that may hinder attempts to accurately characterise such diffusive systems. These errors cannot 
be corrected by increasing the amount of averaged particle locations (either by longer measurements or 
more particles). Rather, they are inherent to noisy, heterogeneous systems. Such systems are the common 
situation in biophysics. As discussed in the introduction, such errors can lead to erroneous conclusions 
about the physical and biological nature of the system. 

Every experimental measurement has a finite accuracy and incorporates some degree of uncertainty. 
In SPT experiments, the location of the tracked particle is known up to a random noise term. As we 
have shown, this added noise causes a systematic decrease of the measured anomalous exponent, even 
when the true MSD is much larger than the noise. For example an MSD to noise ratio of 4 will give an 
error of 20% in estimation of a (equation |4|. We have shown how to quantify and correct this error with 
simple control experiments. 

Heterogeneity in biology may stem from a variety of origins. At the top most level, there is variation 
between individual cells. Even inside each cell, each tracked particle may experience different surround- 
ings, especially in finite time measurements. In addition, biological tracers may have differences in their 
properties, be they artificially inserted, artificially expressed or naturally occurring tracers. 

Many times, the objective of the experimentalist is to characterize the distribution of parameters in a 
system. This is also true regarding the anomalous exponents of the individual particles. However, when 
analysing short tracks, it is difficult to accurately identify the parameters of single particles individually, 
since the small amount of time points introduces large errors. Hence, one aspires to reach conclusions 
from studying ensemble averaged quantities such as the EA-TAMSD. 



8 



We have shown that the variation of anomalous exponents creates a time dependent systematic error 
in the extraction of the average anomalous exponent from the EA-TAMSD. By using the MLSD, we 
are able correct this error and extract the true mean of the distribution. In addition, we study the 
propagation of the errors in the MSD, and manage to quantify the variance of the whole distribution of 
anomalous exponents. 

Finally, we have shown how to combine the noise and the ensemble heterogeneity corrections. Thus 
an accurate MSD of the average particle behaviour is achieved (i.e. the average diffusion coefficient and 
anomalous exponent). This corrected average MSD suffers less random errors due to the large amount 
of points averaged upon and does not suffer the systematic errors of the original EA-TAMSD. Since 
this technique gives both the average anomalous exponent and its variance, the whole distribution of 
anomalous exponents is characterized. 

This ability to characterise the whole distribution from ensemble averaged quantities was shown 
to work even for extremely short trajectories, where single particle analysis failed completely. The 
intuitive way to improve single particle characterization is to measure longer tracks. This is usually 
a complicated experimental challenge which demands more complex and expensive equipment. The 
MSD error estimation technique enables the experimentalist to improve the characterisation of the single 
particle distribution by measuring more particles and not longer trajectories. Thus it carries high potential 
to reduce experimental costs and complexity. 

The deviation of the ensemble averaged NC-MSD (i.e. without heterogeneity correction) from the 
heterogeneity corrected MSD, can give important insight regarding the function of the system. Since the 
ensemble includes particles with high anomalous exponents in addition to particles with low exponents, 
the ensemble covers a larger area than would be expected if all particles behaved according to the average 
particle. At the same time, particles with low anomalous exponents tend to stay more localized around 
the origin of the motion. Thus the first group acts as searchers for distant targets while the second group 
raises the efhciency of interaction near the origin. We see therefore that a distribution of anomalous 
exponents can be beneficial in various sceneries of biological signalling and interactions. 

In the near future, we expect biophysics labs will produce ever growing quantities of data. New tracer 
families, faster acquisition rates and larger quantities of tracers per experiment will enable the develop- 
ment of novel biophysical theories. The correct estimation of parameters from the growing quantities of 
data is essential for these theories to be accurate and preform as useful predictors of biological processes. 

We hope that the techniques described will enable experimentalists to better characterise SPT results. 
This in turn will now doubt improve our understanding of biophysical systems, and help in creating better 
theoretical models for their description. 



Materials and Methods 
Simulation of FBM 

Particle tracks where simulated using the MATLAB wfbm function. Each path i was built for N time 
points with a random Hurst index, Hi, taken from a distribution with mean and variance (for 
FBM, 2H = a). The paths where then down-sampled by a factor of 16 to give T time points. This down- 
samplingremoves high frequency errors in the data. Wfbm uses the algorithm developed by Sellan and 



Meyer 33 
In order 



and thus the simulated diffusion constant of each particle obeys D^ i = T{1 — a) ^"^^^y^^^'* ■ 
'o simulate paths with a minimal correlation between ai and Daj each path was multiplied 
by y/T{l-a)Sin{Tra/2). Since r(l - a)r(l - a) SM^<^mcos{^a/2) ^ ^ ^^^^ significantly reduces the 
correlation. Measurement noise was simulated by adding a random noise increment N^j^f-^ to each track 
at each time point. This increment was taken from a stationary normal distribution with a mean of zero 
and a variance of according to the requested noise scenario. 
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Path analysis 

Since the purpose of this article is to present new methods for the analysis of anomalous processes, 
the analysis of simulated data was performed as if no prior data was available. Specifically, a separate 
Mathematica program was written, which analysed the "experimental" paths only with the help of another 
set of simulated "experimental noise" paths with the same variance and no FBM process. 
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Figure 1: Noise causes a systematic error in TAMSDs. Trajectories of 2^^ time points where 

simulated with a = 0.5 together with an experimental measurement noise of p = 1.2. a) The TAMSD 
(squares) and NC-MSD (circles) of one of these trajectories are shown together with the theoretical 
TAMSD without measurement noise (dashed rising line). The noise TAMSD is shown as a parallel 
continues line. As predicted, the regular TAMSD shows an increasing error at small A while the NC- 
MSD presents a good fit. b)The dynamic functional of the TAMSD (squares) was calculated at various 
time points and was found to follow the predicted erroneous values. The dynamic functional of the 
NC-MSD (circles) on the other hand closely follows the simulated anomalous exponent (continuous line, 
dotted lines are ±10% error of a). Since a finite amount of time points is used for each TAMSD, values 
of vary slightly between individual MSDs. The values presented are averages and perpendicular lines 
show standard errors. 
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Figure 2: Correction of the ensemble variation error, a) For 2^*^ particles, an anomalous exponent 
was taken from a normal distribution with — 0.5 and ctq = 0.2. Total time was T — 2^^ and 2^^ time 
points where taken. b)The TAMSDs of two representative particles are shown. The roughness of the 
path and random deviation from a straight line can be seen, c) The regular EA-TAMSD (circles, blue), 
the MLSD (squares, red) and the theoretical prediction for the average particle (dashed red line) can be 
seen. Lines in MSD and MLSD are to guide the eye. The deviating slope of the MSD can be seen while 
the MLSD is parallel to the theoretical average particle, d) The dynamic functional was calculated for 
the EA-TAMSD (squares, green) and MLSD (circles, blue) at constant intervals on the logarithmic scale. 
The EA-TAMSD dynamic functional is linear and rising (dashed green line), while the MLSD is around 
the simulated fia = 0.5 (blue line). Dotted black lines denote 10% errors from fia- 
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Figure 3: Finding the standard deviation of anomalous exponents. 5000 short trajectories (2^ 
time points, down sampled from paths of length 2^'^), where simulated with fia — 0.6 and CTq, = 0.15. In 
addition an experimental noise of p ^ 1.5 was addded. a) Trying to find the anomalous exponent for 
each track individually fails, as a much wider distribution is found (green) with an ofi^set towards lower 
values, compared to the simulated values (blue), b) For 12 time points, the dynamic functional of the 
EA-TAMSD (square) and MLSD (circle) is calculated. Even though they are both under the simulated 
average (straight line), the differences can be fitted to a straight line on a logarithmic scale (inset). The 
resulting standard deviation found was CTq = 0.16, an error of 7%. 
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Figure 4: Constructing the average particle MSD. 1000 trajectories of 2^" time points where 
simulated from a distribution of anomalous exponents with fi^ — 0.6 and CTq = 0.2. A measurement 
noise of p = 1.5 was added, a) The EA-TAMSD (blue circles) shows strong deviations both in short and 
long times. The corrected average MSD {I'c, black squares) is almost identical to the theoretical average 
particle (dashed red), b) Calculating the dynamic functional for the EA-TAMSD (red squares) and Vc 
(blue circles) shows the dramatic improvement. Blue line is the ensemble average, fia = 0.6, dashed red 
line is the noise induced error and black dotted lines are 10% errors from fia 



